Ultracold polarized Fermi gas at intermediate temperatures 



O 

o 

< 

OO 
(N 



I 

o 
o 



> 

m 
cn 
in 
o 
\o 
o 



X3 

o 
o 



X 



Jani-Petri MartikainerQ 

Department of Physical Sciences, University of Helsinki, 
PO Box 64, 00014 University of Helsinki, Finland 
(Dated: February 6, 2008) 

We consider non-zero temperature properties of the polarized two-component Fermi gas. We point 
out that stable polarized paired states which are more stable than their phase separated counterparts 
with unpolarized superfluid region can exist below the critical temperature. We also solve the system 
behavior in a trap using the local density approximation and find gradually increasing polarization 
in the center of the system as the temperature is increased. However, in the strongly interacting 
region the central polarization increases most rapidly close to the mean-field critical temperature, 
which is known to be substantially higher than the critical temperature for superfluidity. This 
indicates that most of the phase separation occurs in the fluctuation region prior to superfluidity 
and that the polarization in the actual superfluid is modest. 

PACS numbers: 03.75.-b, 32.80.Pj, 03.65.-w 



I. INTRODUCTION 

Since their experimental realization the properties of 
the strongly interacting superfluid Fermi gases 0, 0, S 
H. Isl. 1^. I^.M I9I lioj have been studied ever more actively. 
In particular the recent experiments [ill [T3. [T3 . 
in a polarized two-component Fermi system of trapped 
ultra cold gases have probed the properties of fermion 
pairing in a novel settings. Some important theoretical 
contributions 0, 0| appeared before the breakthrough 
experiments, but the experiments motivated a series of 
theoretical works [H [ll 113, ISl HI HI 111 which have 
explained various features of the experiments with non- 
rotating gases in the low temperature regime. 

Since the BCS-phase cannot support polarization at 
zero temperature, one typically expects, in a harmonic 
trap, phase separation into an unpolarized superfluid 
core, surrounded by the polarized normal Fermi gases. 
The transition between the phases is sharp in the local 
density approximation (LDA) and local minimization of 
the grand potential rules out the possibility of the local 
maxima of the free energy, known as the Sarma 25] or 
breach-pair (BP) [2^ state. LDA approach can, however, 
be modified to include a possibility for a BP phase if one 
minimizes the total (global) energy of the system [2^ . 
Such modification can in principle predict a ring of BP 
phase in the boundary region where BCS phase is trans- 
formed into the normal phase, and for large polarization, 
BP phase is also found in the center of the trap. In 
this paper we will assume that these boundary effects 
are small and can be ignored. For a system in which 
the number of fermions is very large and the trap does 
not have a too high aspect ratio, this is expected to be a 
good approximation although it is quite likely that these 
boundary terms are important in defining the largest pos- 
sible polarization that can still support a superfluid. 
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While the above described works consider mainly zero 
temperature, we now concentrate on finite temperature 
effects. In a work fairly closely related to ours Machida 
et al. considered the finite temperature effects by 
solving the mean-field Bogoliubov-de Gennes equations 
numerically in a cylinder with a harmonic trap in the ra- 
dial direction. Qualitative picture they find is in many 
places similar to ours, although the simple local den- 
sity approximation we apply in a trap is insufficient to 
properly describe the regions they call "FFLO" phase. 
Also a very recent work by Chien et al. "2^ investigated 
the superfluid properties of the homogeneous polarized 
Fermi gas across the BEC-BCS cross-over also at non- 
zero temperatures. At intermediate temperature regime 
they found that stable polarized Fermi superfluids are 
possible. Various aspects of finite temperature effects in 
BEC-BCS cross-over were also recently discussed by Yi 
and Duan [30l| . 

The purpose of this paper is to point out that polar- 
ized BCS states, when they exist, also have a lower energy 
than the corresponding phase separated state with an (fi- 
nite temperature) unpolarized paired part and polarized 
normal part. Furthermore, using LDA we find the state 
of the system in a trap as a function of temperature and 
find that the system becomes gradually polarized in the 
center of the trap. Since this occurs also at temperatures 
still well below the critical temperature we expect the 
mean-field theory employed here to provide at least qual- 
itatively valid description of the system. However, most 
rapid increase of the central polarization occurs in the 
fiuctuation region above the superfluid transition tem- 
perature. This indicates that most of the phase separa- 
tion in the trapped system would occur prior to the onset 
of superfluidity. 

This paper is organized as follows. In Sec. ^1 we dis- 
cuss the uniform BCS phase and compute how highly 
polarized such a state could be. In Sec. lIIII we proceed to 
compare the energies of the polarized BCS phase with the 
mixed phase of unpolarized BCS phase and a polarized 
normal gas. In Sec. II VI we investigate the intermediate 
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temperature properties in a trapped gas using the local 
density approximation. We end with some concluding 
remarks in Sec. 

II. POLARIZED UNIFORM DENSITY BCS 
STATE 

The grand potential corresponding to the standard 
mean-field BCS theory of interacting two component 
Fermi can be computed from the mean-field Hamiltonian 

H = ^[ek-A^.]V3I,A. (1) 

k,(T 

+ ^ H (^L^lk,l + V'k,l'/'-k,2) , 

k ^ 

where is the (fermionic) annihilation operator for 
species a and the coupling strength g — inh^a/m is ex- 
pressed in terms of the scattering length a. Further- 
more, A = g{ijji{Y)ijj2{Y)) is the order parameter/gap 
function, which we take to be independent of position, 
tk = h^k^/2m is the free dispersion, and ii^ is the chem- 
ical potential for the a-component. Since we are going to 
apply this Hamiltonian also in the strongly interacting 
regime, some phenomenological model building is clearly 
involved. In particular, the existence of a non-zero gap 
function should not be taken as an indicator of superflu- 
idity. Superfluidity would imply phase stiffness and this 
typically occurs at a lower temperature Is^l than the 
mean-field critical temperature predicted by the above 
Hamiltonian. 

The assumption of a constant gap function is in this 
context quite reasonable, but it does rule out the possibil- 
ity of the FFLO-phases [33, |3J| , which could in principle 
exist in some narrow range of parameters ■ Diagonal- 
izing the above Hamiltonian we can simply calculate the 
grand canonical free energy 

n{^il,^l2,A) = -fcBriog[rr(exp(-/3i7))], (2) 

as a function of temperature T |23| • Due to the use of the 
contact interaction the resulting expression is ultraviolet 
divergent and this divergence is removed by the usual 
renormalization 

A2 A2 

where V is the quantization volume which drops out 
from the final physical results. The gap equation can 
be derived from the grand potential as condition for the 
extrema dCl/dA = 0, while the number equation ap- 
pears from rio- = —dQ/dfj,a, where is the component 
density. Alternatively, the latter condition can be en- 
forced as an extremal point of the Helmholtz free energy 
F = Q + fiiUi + ^2^2, i-e. dF/dficr = 0. In this section 
we choose the units using the ideal Fermi gas of density 
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FIG. 1: The highest value of polarization that the uniform 
BCS state can support as a function of temperature (in units 
of Fermi temperature) when kpa — —100 (solid line) and 
kpa = —1 (dashed line). We have set the maximum polariza- 
tion equal to zero in the region above the critical temperature. 



n as a benchmark and h is taken to be an average of the 
component densities. This means that the unit of length 
is given by Z = l/kp — {6n^n)^^/'^ , while the unit of 
energy is ep = h^kp/2m. 

The polarized BCS solution corresponds to the min- 
imum of the grand canonical free energy, but does not 
exist at zero temperature. At zero temperature only sta- 
tionary points of the free energy corresponding to the 
polarized solutions and spatially constant gap are the 
normal state and the BP state. However, as the tem- 
perature is increased stable polarized BCS solutions can 
appear. 

In Fig. n we show the largest possible polarization 
Pc{T) — {ni — n2)/(ni -|- 712) with a large and moder- 
ate coupling strength. It is clear that BCS solutions cor- 
responding to the minimum of the grand potential with 
substantial polarizations do exist in the strongly interact- 
ing regime as the temperature increases. For very weak 
interactions the polarization supported becomes vanish- 
ingly small, but the maximum polarization of the BCS 
state can be as high as Pc{T) « 0.68 in the strongly 
interacting regime with kpa — —100. If the polariza- 
tion is higher than this the system will never enter a 
uniform superfluid region as the temperature is lowered 
towards T — 0. Curiously, this critical density asymme- 
try is close to the experiment ally observed critical atom 
number asymmetry in a trap |l5| . 

Naturally close to critical temperature of the mean- 
field theory the approach used here becomes unreliable, 
but states of considerable polarization do exist even well 
below the critical temperature. The maximum polariza- 
tion supported by the BCS state is predicted to peak in 
the intermediate temperature regime. This is expected 
since the critical temperature is maximized in an unpo- 
larized system and therefore just below the critical tem- 
perature only small polarization is possible. 
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III. PHASE SEPARATED STATE VS. 
POLARIZED BCS STATE 

Let us now turn to the possibility of phase separation. 
At zero temperature the phase separated state is clearly 
favored since the polarized BCS solution does not exist. 
However, as we saw in the previous section at non-zero 
temperature the polarized BCS solutions do exist and it 
is not clear if these polarized BCS states have a lower 
energy than the phase separated state corresponding to 
the same number of atoms. 

In particular, we have in mind a comparison with an 
unpolarized BCS phase occupying a volume fraction 1 — 
X {0 < X < 1) next to a polarized normal phase in a 
volume fraction x 'sF]. If the unpolarized BCS phase 
has a density nscs the normal state will have fermion 
densities 



- (1 - x)nBcs 



(4) 
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FIG. 2: The optimal volume fraction x for the mixed state 
of the unpolarized BCS phase and normal polarized gas at 
T = as a function of polarization for kpa = —1 (dashed) 
kpa — —10 (solid). 



if the average density for each component, averaged over 
the whole the system, is Ua-. The problem is then 
to find the values of x and ubcs which minimize the 
Helmholtz free energy F{x,nBcs) for given set of pa- 
rameters, such as coupling strength, polarization, and 
temperature. When the interface energy is ignored the 
Helmholtz free energy is simply a sum of two contri- 
butions, the polarized normal gas and unpolarized BCS 
phase. 

In Fig. 121 we show how the optimum normal state vol- 
ume fraction behaves as a function of polarization in a 
phase separated state at zero temperature. Quite sen- 
sibly, for weaker interactions the normal part increases 
faster with polarization while in the strongly interacting 
regime the paired part persists longer. Our result for 
weaker interaction appears to be consistent with the ear- 
lier result by Bedaque et al. 35]. As we increase the tem- 
perature the normal part will become bigger and eventu- 
ally at the critical temperature we have a phase transition 
into a pure normal state. 

Since our approach does not include surface energy 
contribution from the boundary between superfluid and 
normal regions, the zero temperature transition to a 
fully normal state as polarization increases is not sharp. 
Nevertheless, the volume fraction of the paired region 
becomes so close to zero that it drops below our nu- 
merical resolution. If one were to include the surface 
terms Is^ , one would expect a sharp boundary as well as 
a shift downwards from the "critical" polarization around 
p w 0.85 for the stronger interactions. 

In Fig. 121 we show the minimum attainable Helmholtz 
free energy as a function of temperature for the phase sep- 
arated state, together with the energies of the polarized 
BCS state, and the pure normal state with fixed polariza- 
tion of 0.1 in the strongly interacting regime. While the 
energy differences are not large, it is nevertheless clear 
that whenever the (stable) polarized paired state exists, 
it is more stable than the phase separated state. Fur- 




FIG. 3: The free energy of the uniform polarized state (solid 
line), phase separated state (dashed line), and the polarized 
normal state (dash-dotted line) as a function of temperature 
(in units of Fermi temperature). We choose the polarization 
p = 0.1 and coupling kpa = —20. The line for the uniform po- 
larized state begins from the temperature where the uniform 
polarized phase is first possible. 



thermore, the energy of the phase separated state is an 
underestimate due to the absence of the surface energy 
contributions. The mean-field theory used here ignores 
some fluctuations contributions to the free energy, but we 
believe these are not going to change the outcome of the 
energy comparison since these contributions would have 
to be included in both competing states. This outcome is 
made all the more likely by the absence of surface energy 
contributions in the phase separated state. 
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IV. LOCAL DENSITY APPROXIMATION AT 

T / 

At zero temperature, for an ideal Fermi gas we find in 
the local density approximation a component density 



rir) 



1 

6^ 



2™ 



-1 3/2 



Vrir)) 



(5) 



where the external trapping potential is given by Vri'r) — 
mup'r^ 12 and the component chemical potential 
is related to the particle number N„ through /i^ = 
'hu:{%N„Yl'^ . The density vanishes outside the Thomas- 
Fermi radius RTF,a = \/^4i^aJ™j^ ■ The polarization 
Vn — in\ — n2)l{n\ + ^2) will then be 

(a^i - VT(r))3/2 - (^2 - yT(r))3/2 



VNif) 



(m-^T(r))3/2 + (^2-^^T(r))3/2 



(6) 



inside the Thomas-Fermi radius of the minority compo- 
nent. This implies that for an ideal gas the polarization 
will increase from 



PAr(r = 0) 



1 + V^V^Vi 



(7) 



in the center of the trap to the maximum value of 
PA'(r > Rtf.'i) = 1 beyond the Thomas-Fermi radius 
of the minority component. 

On the other hand at zero temperature superfluid one 
expects strong phase separation due to the absence of the 
polarized superfluid state corresponding to the minimum 
of the grand potential. In that case, if we ignore the 
boundary effects, we would find vanishing polarization 
Ps(r) in the center of the cloud, which then abruptly 
increases close to the trap edge, eventually becoming 
Psir > Rtf,2) = 1 beyond the Thomas-Fermi radius 
of the minority component |23j. 

As we have found, the polarized paired state can exist 
at non-zero temperature and furthermore when it exists 
it has a lower energy than the corresponding phase sep- 
arated state with non-polarized superfluid mixed with a 
normal Fermi gas. This then raises the question of what 
the polarization p{r) will be in a trap at intermediate 
temperatures. Finite temperature effects will round off 
the density distribution and make the LDA analysis less 
accurate, but we expect it to be qualitatively valid. Also, 
since the mean-field theory of the strongly interacting 
Fermi gas becomes more unreliable at higher tempera- 
tures, it appears fruitless to aim towards higher precision 
within the framework of the mean-field theory. 

Since the energy of the polarized BCS state increases 
with polarization it will be energetically more favorable 
to have a large number of atoms in a region of fairly 
small polarization, rather than forming, for example, a 
small core of highly polarized superfiuid surrounded by 
the normal gas or perhaps unpolarized BCS phase. This 
implies, quite sensibly, that p{r) wll be something in be- 
tween the normal state polarization pjv('') and the T = 
phase separated superfluid polarization ps(r). 



The normal state polarization increases smoothly up 
to 1 and as we have found, the BCS state can be polar- 
ized if the (local) polarization is below the critical value 
Pc{T). Also one expects that the extent of the (possibly) 
polarized superfiuid core becomes smaller as the temper- 
ature is increased and eventually disappears at a critical 
temperature where polarized BCS states can no longer 
be found. Outside the paired core one expects a polar- 
ized normal Fermi gas. These qualitative arguments are 
quite well confirmed by solving the state of the system 
using LDA analysis. This implies using the local chemi- 
cal potentials /ie//,o- = Mo- — Vrir) and finding /io- in such 
a way that the total particle numbers have the desired 
values. At non-zero temperatures the contributions from 
/^e//,CT < do not vanish identically and the edge of the 
cloud is rounded somewhat by this effect. 

In Fig. 0] we show a typical result for the densities, the 
gap, as well as the density difference as a function of po- 
sition and temperature in the strongly interacting regime 
with a number asymmetry (A^i — N2)/{Ni + N2) = 0.5. 
The behavior of the gap function is qualitatively simi- 
lar even if the paired region were polarized. The spatial 
extent of the paired core first, somewhat surprisingly, in- 
creases slightly as the temperature increases, but is then 
reduced as one approaches the critical temperature. This 
effect is however so small, that it would be difficult to 
observe and could be an artifact of the LDA. At zero 
temperature LDA analysis gives rise to abrupt changes 
in densities (especially for the minority component) and 
also the gap drops suddenly to zero. Increasing tem- 
perature eventually removes all these discontinuities. As 
the temperature is lowered below the critical tempera- 
ture the minority component develops a clear bimodal 
density distribution. 

The polarization as a function of position naturally be- 
comes different as the temperature increases due to the 
appearance of polarization in the paired core. The po- 
larization in the center of the cloud increases gradually 
with increasing temperature, but this build-up becomes 
more abrupt if the number asymmetry increases. If fiuc- 
tuation effects were to be included, the critical temper- 
ature for the superfluidity could become so low (in 
the strongly interaction region) that the actual superfiuid 
region would have only a small polarization. The build- 
up of polarization and removal of phase separation with 
increasing temperature would therefore mostly occur in 
the "pseudo-gap" region. For more weakly interacting 
system this distinction disappears, but then also the po- 
larization that the BCS state can carry drops and the 
superfiuid region will again have only a small polariza- 
tion. 

Qualitatively the behavior predicted by the LDA 
does not change strongly when the number asymme- 
try is increased to 0.8. Interestingly, it turns out 
that the position dependence of the local polarization 
[ni{r) ~ n2{2)] / [ni{r) + n2{r)] becomes weaker as the 
temperature increases. This suggests that approximat- 
ing the system as uniform, might not always be too dras- 
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FIG. 4: (a) The pairing gap, (b)-(c) component densities, and 
(d) the density difference as a function of position and tem- 
perature. The result is calculated in the strongly interacting 
regime with coupling strength kpa = —20 and number asym- 
metry (A^i - A^2)/(A''i + A^2) = 0.5. The unit of energy is 
the Fermi energy ep = fLkp/2m of the ideal Fermi gas of ^Li 
atoms with (A'^i + N2)/2 = 10^ atoms in a harmonic trap with 
2n{ujxUjyUJz)^^^ = 600 Hz. Unit of length is l/kp and densities 
are scaled with the factor Gtv^ so that ideal gas would have a 
density of riideai = 1 in the center. 



tic approximation, especially in the neighborhood of the 
temperature when pairing effects begin. 

Since it is experimentally easier to measure inte- 
grated densities such as the doubly integrated na-{z) — 
J dxdyn„{x,y, z) or the column integrated na-{z,y = 
0) = / dxncr(x, y — 0, z) we show in Fig.[Slthe integrated 
density differences for the polarized and unpolarized su- 
perfluid systems at zero as well as non-zero temperature. 
As can be seen the signals are not qualitatively very dif- 
ferent. The polarized state has somewhat higher central 
density difference and finite temperature effects make the 
profiles of the polarized system more smooth. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have employed the mean-field theory 
at finite temperatures to discuss the behavior of the inter- 
acting Fermi gases with unequal numbers of fermions in 
two components. We found that stable polarized paired 
states can exist at intermediate temperatures and that 
cross-over to a polarized paired state can leave a mod- 
erate signature on the integrated density distribution of 
the Fermi gas. The LDA approach used here predicts 
a gradual appearance of the polarized superfiuid as the 



temperature is increased. 

The accuracy of the mean-field theory in the strongly 
interacting regime is naturally questionable especially at 
non-zero temperatures and is expected to give rise, for 




FIG. 5: (a) The doubly integrated density difference An{z) = 
/ dxdy [ni(r) - 712 (r)] at T = (solid) and T = 0.2 (dashed), 
(b) The column integrated density difference An(z, y — 0) = 
J dx [ni{x,y = 0, z) — n2{x,y = 0, z)] at T = (solid) and 
T — 0.2 (dashed). The result is calculated with same pa- 
rameters and units as in Fig. 3] 



example, to a substantial overestimate of the critical tem- 
perature for superfiuidity Therefore, the critical 
temperature found here provides an estimate on the tem- 
perature when pairing effects appear, but not necessar- 
ily on the superfiuid transition temperature when phase 
coherence sets in [s^. The results presented here in the 
strongly interacting region are qualitatively similar in the 
weakly interacting regime where the fluctuations are less 
strong. The absolute magnitudes of the possible polar- 
izations, the superfluid gap, and the critical temperature 
are then, however, substantially smaller. 

At zero temperature it is known from the Monte-Carlo 
calculations that the mean-field theory overestimates 
the gap substantially. In order to provide more accurate 
quantitative predictions it would therefore be important 
to understand the intermediate temperature properties 
of strongly interacting Fermi gases better. This is all 
the more important due to the very recent MIT exper- 
iment |l5j | which observed prominent interaction effects 
close to and somewhat above the (superfiuid) transition 
temperature and also pointed out the need of incorpo- 
rating the interaction effects between the normal and su- 
perfluid regions into a quantitatively accurate theory. 
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